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Is there an A y problem in low-energy neutron-proton scattering? 
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We calculate A y in neutron-proton scattering for the interactions models WJC-1 and WJC-2 
in the Covariant Spectator Theory. We find that the recent 12 MeV measurements performed at 
TUNL are in better agreement with our results than with the Nijmegen Phase Shift Analysis of 
1993, and after reviewing the low-energy data, conclude that there is no A y problem in low-energy 
np scattering. 



In a recent paper, Braun et al. (designated BR08) re- 
port on a new set of very precise measurements of A y 
in elastic neutron-proton scattering at 12.0 MeV neutron 
lab energy |T] . They compare their experimental observ- 
ables to the results produced by the Nijmegen Phase- 
Shift Analysis [2 of 1993 (PWA93) and by the CD-Bonn 
potential model [3], also fitted to the data and producing 
phase shifts consistent with PWA93. They observe that 
the predictions of the PWA93 as well as the CD-Bonn dif- 
fer considerably from the new data and interpret this as 
evidence for a significant shortcoming of the phase shifts 
and, consequently, of all potential models fitted to them. 
In a model study, starting from the CD-Bonn potential, 
they find that increasing the charged-pion coupling con- 
stant 5^ ± /47r, while keeping the neutral-pion coupling- 
constant unchanged at g 2 /47r = 13.6, brings the the- 
oretical predictions into better agreement with the new 
data. 

In this letter, we take up the issues raised in [T] and 
consider them in light of a recently finished new analy- 
sis of np scattering within the Covariant Spectator The- 
ory (CST)H 0. As described in detail in Refs. [6], 
the scattering amplitude M for two-nucleon scattering in 
CST is obtained by solving a covariant integral equation 
of the form 

M = V — VGM (1) 

where V is the irreducible kernel (playing the role of a po- 
tential) and G is the intermediate state propagator. The 
Bethe-Salpeter (BS) equation [7j has a similar structure, 
but in the BS theory the four-momenta of the particles in 
intermediate states are subject only to the conservation 
of total four- momentum P = P1+P2, so loop integrations 
over intermediate states are over four independent vari- 
ables. The CST equation maintains four-momentum con- 
servation, but constrains the energy components by plac- 
ing one of the two intermediate particles on its positive- 
energy mass shell. Since P is fixed, both intermediate en- 
ergies become functions of the three-momenta only, and 



TABLE I: Statistical errors for the data set BR08 of Ref. [J. The 
systematic error is 1.5%. 
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all intermediate loop integrations reduce to three dimen- 
sions, as in the non-relativistic theory In spite of this 
considerable simplification, the equation remains covari- 
ant because the on-mass-shell constraint itself is a covari- 
ant condition. This framework has been generalized to 
other systems and applied successfully to many problems, 
in particular also to the three-nucleon system |8j [9] . 

Our recent fit to the np scattering data [5] uses the 
CST with a kernel approximated by the sum of one- 
boson-exchange (OBE) interactions. Our best model, 
designated WJC-1, uses 8 bosons with 27 parameters 
and provides a high precision fit to the 2007 data base 
(x 2 /iVdata = 1.06). A new phase shift analysis emerges 
from this fit, with phase shifts that differ significantly in 
some cases from the PWA93 analysis. Another model 
(WJC-2) was chosen to be as simple as possible, and re- 
quires only 6 bosons with 15 parameters. It does almost 
as well, with xV^data = 1.12. 

Before making a detailed assessment, we point out that 
the errors quoted by BR08 combined systematic and sta- 
tistical errors together in quadrature. A better procedure 
is to treat the systematic normalization error as an inde- 
pendent experimental degree of freedom, so that the \ 2 
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FIG. 1: Neutron-proton analyzing power at 12.0 MeV. Lines are WJC-1 (solid), WJC-2 (dashed) and Nijmegen PWA93 (dotted). Left 
panel shows BR08 data scaled by WJC-1 (solid circles) and WJC-2 (half filled squares). Right panel shows the WE92 |10| data scaled by 
WJC-1 (half filled squares) and BR08 data scaled by PWA93 (solid circles), 



for the data set would be 
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where Oi and t, are the measured and the calculated value 
of the observable at point i, 5oi and <5 sys are the statistical 
errors at point i and the systematic error, Z is a factor 
[chosen to minimize Eq. ([2])] by which the data and errors 
can be divided to correct for the systematic error and 
improve the agreement with theory, and the last term is 
the additional contribution to the x 2 coming from the 
renormalization of the data. We contacted the authors 
of BR08 who told us that the systematic error was about 
1.5%, and provided us with the original statistical errors, 
reproduced in Table [I] [FT] . 

Using the new errors given in Table [T] and the new 
models (and phase shift analyses), we reach a different 
assessment of the impact of the BR08 experiment. Fig- 
ures [T] and [2] present different comparisons of the 12 MeV 
A y BR08 measurements with the predictions of the Ni- 
jmegen PWA93, and Models WJC-1 and WJC-2. Table 
iT| gives the X 2 /-^data for each of these theories and Table 



III surveys the previous A y measurements for all labora- 



tory energies less than 20 MeV. From these figures and 
tables we draw the following conclusions: 

First, without even looking at the data, observe that 
the predictions of both of our models WJC-1 and WJC-2 
differ substantially from the PWA93 prediction (the the- 
oretical curves shown in both panels of Fig. [IJ are iden- 
tical), expecially at the smaller angles. This is not un- 
expected; it is a consequence of the difference between 
our phase shifts and those of PWA93. This may have 



TABLE II: The average \ 2 from Ref. pQ compared to the average 
X obtained by scaling the data and errors in Table [I] [with the 
scaling factor Z determined by minimizing Eq. |2j]. The scaled 
WE92 I10| data, also measured at 12 MeV, are shown for compari- 
son. The Nijmegen rejection criteria say that a set with this many 
data (15 points plus one normalization error) should be retained 
only if x 2 /-^data < 2.26. Note that the BR08 data almost meets 
this criteria with model WJC-2. 



Model 
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BR-scaled 


Z(BR) 


WE92 


PWA93(l) a 


6.621 


3.849 


0.9367 


1.71 


PWA93(2) a 


6.654 


3.861 


0.9365 


1.71 


WJC-1 


3.631 


3.014 


0.9690 


1.61 


WJC-2 


2.377 


2.387 


1.0029 


1.16 



"Our fitting procedure uses the effective range expansion. The 
numbers shown for PWA93(2) use WJC-1 parameters, which give 
the best fit to the overall data set at all energies. Nijmegen does 
not give 1 Sq parameters, but PWA93(1) uses 3 Si parameters taken 
from Ref. 1121. which does not fit the overall data set as well. 



significant impact on other observables as well. 

Next, Table [TT] shows that the scale factor Z differs 
for each model, and hence there is no one "correct" data 
set. For example, to accommodate PWA93, the data is 
increased by over 6% (four standard deviations), giving a 
large normalization contribution to the x 2 but lowering 
the collective effect of the statistical errors on each of the 
15 measured points. This renormalization of the data 
lowers the X 2 /A 7 d a ta from about 6.6 to about 3.85. In the 
end, even with this renormalization the PWA93 model is 
the least successful in describing the data. For compar- 
ison, Fig. [TJ (right panel) and Table [TT] also show previ- 
ous A y measurements performed at TUNL at the same 
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FIG. 2: Ratios of the 12.0 MeV BR08 data for the analyzing power 
to the predictions of WJC-1 (1), WJC-2 (2), and PWA93 (N). 



energy, designated WE92 [ID]. Qualitatively, they ap- 
pear to agree well with the new measurements, and also 
with the theoretical models, with considerably smaller 
X 2 /iV<iata than BR08, owing mainly to their larger statis- 
tical errors. 

How are we to interpret these large BR08 x 2 ? Does it 
mean that all of the theories are wrong? In this context 
we must remember that a phase shift analysis gives the 
most general theoretical description of scattering below 
the pion production threshold, and while several theoret- 
ical models might generate the same set of phase shifts, 
the best phase shifts provide the best model independent 
description of the data. At low energies, only a few par- 
tial waves contribute noticeably to the NN scattering am- 
plitude, and of those only the S-waves are not small. In 
principle, this leaves the remaining partial waves with 
large uncertainties. However, the low-energy phase shifts 
arc linked to the much better determined phases at higher 
energies through the requirement of a smooth continuous 
energy dependence, which limits their uncertainties and 
relates observables at different energies. While a contin- 
uous energy dependence arises automatically when a po- 
tential model is used, it is imposed also in a multi-energy 
phase shift analysis, although in a phcnomcnological way. 

Thus, any new data set effectively competes with the 
entire phase shift analysis, which is based in turn on all of 
the 3788 data currently included in the 2007 data base. 
If a new data set disagrees with the predictions of the 
phase shift analysis it also disagrees with all of the other 
data. 

So the BR08 data disagrees with the latest phase shift 
analyses (and therefore the entire 2007 np data base), 
but does this represent a real incompatibility? Is the 
xV-^Vdata of this data set too high to be tolerated in a 



quantifiable sense? The Nijmegen group used a statisti- 
cal criteria based on the following consideration: from a 
given model (in our case a phase shift analysis or a po- 
tential model) one calculates the observables contained 
in a measured dataset. If the model described the ob- 
servables exactly, in other words, if we had the "correct" 
theory, one should obtain x 2 /-^data ~ 1 as a result of the 
unavoidable and random statistical errors in the experi- 
mental data. This presupposes that each individual data 
point in the set is normally distributed around the theo- 
retically predicted value with a standard deviation equal 
to the statistical error. 

One can then ask: based on the theoretical model, 
and allowing for random fluctuations of the size given by 
statistical errors, what is the probability that one will 
obtain a xV-^data as large as or larger than the actu- 
ally observed value for the data set? If that probability 
becomes very small, for instance smaller than some pre- 
defined value, one may conclude that the model and the 
data are not compatible with each other. By the same 
line of arguments, getting very small values of X 2 /-^data 
can also be too improbable, leading to a second criteria 
testing data sets for too low xV-^data- 

The Nijmegen criteria adopts a critical probability of 
about 0.27%, in analogy to the usual "3cr-criteria" for 
testing individual normally distributed measurements. 
Accordingly, for A^ data = 16 (15 data points plus one com- 
mon normalization error), the highest admissible value of 
xV^data is 2.26. From Table [n] we see that the BR08 
data set is not compatible, in this sense, with any of the 
models, although it comes very close to WJC-2. 

How significant is this incompatibility, and could it 
be due, in part, to a problem with the BR08 data them- 
selves? Figure 2]shows the ratio of the BR08 A y measure- 
ments to the theory for three different theories. There is 
a rather pronounced disagreement between the data and 
all models at smaller angles, but the "break" in the set at 
about 110° may be more significant. The theories are in 
good agreement with each other at these larger angles, 
yet none of them can reproduce this behavior. Elimi- 
nating, for example, the two data points at 112° and 
120° would already bring the xV^data of BR08 with re- 
spect to both WJC-1 and WJC-2 below the critical value 
(although not for PWA93). A small increase in the es- 
timated statistical errors would produce a similar result. 
It is not our intention here to second-guess the data anal- 
ysis of pQ, but rather to point out that the disagreement 
of their data with the WJC models is mild and may be 
partly due to the behavior of the data set itself. 

In this context it is interesting to examine the np 
database for any other independent evidence of a problem 



with Ay at low energy. Table III shows all A y data for 



E\ a b < 20 MeV. All of these sets (except the backward 
angle T088 measurements which are excluded because 
the x 2 is t°° small) are accepted by the Nijmegen cri- 
teria, and all agree well with both WJC models. Thus 
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TABLE III: Previous A v (or P) measurements for i?i a b < 20 
MeV. The tabulated (x 2 / N )i are for model WJC-1; (x 2 /A r )2 are 
for model WJC-2. 



TABLE IV: Values of the pion coupling constant used in the fits 
described in this paper. 



Slab 


ref 


N 


(x 2 /^v)i 


(X 2 /N)2 


(X 2 /^)PWA93 


7.6 


WE92 


5 


2.10 


2.23 


2.51 


10.0 


H088 


13 


0.76 


1.25 


0.82 


11.0 


MU71 


1 


0.04 


0.13 


0.03 


12.0 


WE92 


9 


1.61 


1.16 


1.71 


13.5 


T077 


1 


0.04 


0.53 


0.02 


14.1 


BR81 


11 


0.36 


0.42 


0.37 




WE92 


6 


0.63 


0.66 


0.67 


14.5 


FI77 


9 


0.87 


0.88 


0.99 


14.8 


T077 


1 


0.35 


1.67 


0.41 


16.0 


T077 


1 


0.04 


0.28 


0.03 




WE92 


6 


1.20 


1.39 


1.10 


16.2 


GA72 


3 


0.17 


0.18 


0.18 


16.4 


BE62 


4 


0.69 


0.71 


0.69 




J074 


4 


0.90 


1.08 


0.89 


16.8 


MU71 


1 


0.03 


0.07 


0.02 


16.9 


M074 


5 


0.63 


0.52 


0.59 




T088" 


12 


1.30 


1.30 


1.33 




T088 h 


5 


0.03 


0.05 


0.06 


17.0 


WI84 


7 


0.52 


0.55 


0.57 


18.5 


WE92 


5 


0.56 


0.47 


0.62 


19.0 


WI84 


7 


0.62 


0.62 


0.62 




All 


111 


0.85 


0.90 


0.90 



"This set of forward angle measurements from 51.0 to 143.7 de- 
retained in the fit. 

set of backward angle measurements from 
degrees is excluded from the fit because Us error is 



grces is 
'This 



136.5 to 166.5 
too small. 



there is no indication of a problem in the older A y data. 

One might argue that the previous A y measurements 
were not as precise as the BR08 data, for instance because 
they were not corrected for the polarization dependent 
efficiency of the neutron detectors pQ. Moreover, they 
have larger errors and are therefore not as restrictive as 
the new precise BR08 data. This may be true, but it 
also means that the larger \ 2 of BR08 is not a result of 
the older A y data pulling the fits into a wrong direction, 
since they could accomodate larger variations in the fits 
without a prohibitive increase in the overall \' ''■ The fits 
must be dominated by other data. 

Our np data base contains 73 data sets below 20 MeV. 
Only three of them are excluded after applying the Ni- 
jmegen criteria: two because their \ 2 is too small, and 
only a single set, consisting of three total cross section 
measurements between 0.5 and 2.0 MeV, because of a 
very large \ 2 - So again, we find no hints for any prob- 
lems in the database. On the contrary, the low-energy 
data seem to be remarkably consistent with each other 
and with the phase shift analyses. 

This brings us back to the question if the BR08 data 



coupling 


WJC-1 


WJC-2 


PWA93 


9lo /4vr 


14.608 


14.038 


13.567 




13.703 


14.038 


13.567 



set should be excluded from further fits because of its 
incompatibility with the rest of the data base. From the 
description given above, it should be clear that, although 
the Nijmegen criteria leads to a well defined critical value 
for the maximum tolerated x 2 /-/V<j ata of a data set, its 
choice is also somewhat arbitrary. 

Instead of worrying too much about whether or not 
the Nijmegen criteria should be relaxed a bit in order to 
accomodate the BR08 data set, we decided simply to see 
if its inclusion in the data base leads to any significant 
changes in the resulting fits. We found that, because it 
adds only 16 points (15 plus the normalization error) to a 
data base of over 3700 points, it had essentially no effect 
on the fitting process. It is important to realize that the 
other 3700 points already fix the phases to a very large 
extent. We must accept the fact that the BR08 data set 
has a large x 2 that cannot be reduced by further fitting. 

We are left with some ambiguities. We cannot con- 
clude, as is suggested in [T], that the BR08 data point 
to a need to increase the pion coupling constants. As it 
turns out, both of the WJC models have larger couplings 



than advocated by the Nijmegen group (see Table IV I, 
but the BR08 data were not used to obtain these results, 
and these results are not altered by including the BR08 
data in a refit. Still, this observation may explain why 
the BR08 data is closer to models WJC. 

One cannot escape the fact that it is inappropriate to 
draw strong conclusions from a single data set. As new 
data sets are added to the data base, they may slowly 
change the phase shifts, but any new data set cannot be 
expected to exert much leverage, and if a new data set 
disagrees significantly with the phase shift analyses we 
are rather led to look for problems with the new data 
set. Our conclusion at the moment is that the BR08 
dataset is marginally consistent with the WJC-2 phases, 
less consistent with the WJC-1 phases (the best so far), 
and probably inconsistent with the old PWA93 phases. 

Acknowledgements: We thank W. Tornow for help- 
ful correspondence and for providing us with the errors in 
Table [T] This work is supported by Jefferson Science As- 
sociates, LLC under U.S. DOE Contract No. DE-AC05- 
06OR23177. A. S. was supported by FCT under grant 
No. POCTI/ISFL/2/275. 



[1] R.T. Braun et al., Phys. Lett. B 660 (2008) 161. 
[2] V.G.J. Stoks et al., Phys. Rev. C 48 (1993) 792. 



[3] R. Machleidt, Phys. Rev. C 63 (2001) 024001. 
[4] F. Gross and A. Stadler, Phys. Lett. B 657 (2007) 176. 
[5] F. Gross and A. Stadler, Phys. Rev. C 78 (2008) 014005. 
[6] F. Gross, J.W. Van Orden and K. Holinde, Phys. Rev. 

C 45 (1992) 2094. 
[7] E.E. Salpeter and H.A. Bethe, Phys. Rev. 84 (1951) 1232. 
[8] A. Stadler and F. Gross, Phys. Rev. Lett. 78 (1997) 26. 
[9] A. Stadler, F. Gross and M. Frank, Phys. Rev. C 56 

(1997) 2396. 

[10] G.J. Weisel et al., Phys. Rev. C 46 (1992) 1599. 
[11] W. Tornow, private communication (2008). 
[12] J.J. de Swart, C.P.F. Terheggen and V.G.J. Stoks, 
arXiv:nucl-th/9509032 (1995). 



